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1 Introduction 



The factorization method has been historically atributed to Schrodinger [^. This ap- 
proach (originated from the well known treatment of the harmonic oscillator in non- 
relativistic Quantum Mechanics) avoids the use of cumbersome mathematical tools. Since 
Schrodinger the factorization method has been successfully applied to solve essentially any 
problem for which there exists exact solution 0. Considering the narrow set of exactly 
solvable potentials, it is interesting to look for new potentials for which the corresponding 
Schrodinger equation becomes analitically solvable. There are a lot of results on the mat- 
ter, in particular it is worthwhile mentionning a work of Mielnik , where a variant of the 
standard factorization was introduced in order to get a family of Hamiltonians isospectral 
to the harmonic oscillator. The first application of Mielnik's factorization was performed 
for the hydrogen-like radial potentials the resulting family of Hamiltonians became 
isospectral to the standard radial hydrogen-like Hamiltonians. An additional interesting 
application arose from the arena of the coherent states. Concerning this subject, it is 
well known that the usual first order factorization operators for the harmonic oscillator 
Hamiltonian are also its creation and anihilation operators. This is not the case for the 
isospectral oscillator Hamiltonians derived by Mielnik because their first order factoriza- 
tion operators are different from their creation and anihilation operators, which turn out 
to be of order equal to (or greater than) three ^. These last operators are the genera- 
tors of non-linear algebras, and allow one to get immediately the corresponding coherent 
states as the eigenstates of the annihilation operator [Q. Recently, these algebras have 
been rederived (using different methods) and reinterpreted by other authors [§,0]. It has 
been shown as well that such non-linear algebras can be associated to other systems 0. 

The factorization techniques 0, |[, the Darboux transformation [§, the super- 
symmetric Quantum Mechanics (SUSY QM) j^, and other procedures have a common 
point: they can be embraced in the algebraic scheme known as first order intertwining 
transformation which makes use of a first order differential operator to intertwine (the 
eigenfunctions of) two different Hamiltonians [0. The generalization of this technique 
involves higher order differential operators and leads in a natural way to the higher order 
SUSY QM [11-15]. One of the advantages of the intertwining technique is that the com- 
position of n first order intertwining transformations allows one to express the nth-oider 
intertwiner operators as the product of n first order intertwining operators |ll|, |13|, |15| . 

The experimental study of spectra (for systems as nuclei, atoms, molecules, etc) makes 
important per se the theoretical determination of the energy levels of potentials capable 
of being used as models for describing different physical situations. However, for a quan- 
tum system with a given Hamiltonian it is seldom possible to get exact solutions to the 
corresponding eigenproblem. 

In this paper we are going to show that, given an exactly solvable potential as an 
input, the intertwining technique provides, in general, new exactly solvable potentials 
which, from the traditional viewpoint, would be solved by perturbative or approximate 
methods. Hence, working with Hamiltonians whose spectrum is known, one has at hand 
the tools to test the validity and convergence of some perturbative or approximate meth- 
ods |l^, |T^. In Section 2 we shall present the nth-order intertwining technique, discussing 
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in more detail those aspects which have not been previously reported in the literature. 
In particular, in Section 2.1 we shall show that the standard Infeld-HuU factorization as 
well as the modified Mielnik's method can be recovered by the first order intertwining 
treatment. We will introduce a generalization of the Mielnik's approach by using solu- 
tions to the Schrodinger equation corresponding to factorization energies not belonging 
to the physical spectrum of the initial problem. In Section 2.2, by using the second order 
intertwining approach, we will introduce explicit expressions allowing the construction of 
two-parametric families of potentials whose spectrum is equal to the initial spectrum plus 
two new levels at predetermined positions. In Section 2.3 we shall prove that the second 
order, as well as the higher order cases, can be recovered by the iterative application of 
our first order intertwining procedures. Section 3 contains the derivation of a wide set 
of new ra-parametric families of potentials whose spectrum is the same (isospectral), or 
almost the same (almost isospectral), as the corresponding hydrogen-like spectrum. We 
shall introduce a set of factorization energies generalizing the choice made for the Infeld- 
HuU and Mielnik factorizations of the hydrogen-like potentials. As a consequence, the 
first n energy levels of the ri-parametric families can be placed at any previously fixed 
spectral positions chosen from the above mentionated set of factorization energies. This 
result opens the possibility of placing holes between these levels or between them and the 
other ones. Some particular cases related to results previously derived by other authors 
are also mentioned. 



2 The Intertwining Method 

Let H and H be the one dimensional Hamiltonians 

H^-f-^+V{x), H^-f- + V{x). (1) 

The nth-oTdei intertwining method aims to find a nth-order differential operator A such 
that 

HA = AH. (2) 

Let if be a Hamiltonian with known solutions to the time-independent Schrodinger 
equation Hip = Eip. We are looking for solutions to the corresponding equation for 
the intertwined Hamiltonian H, Hip = Eip. The transformation (|^) preserves many 
spectral properties of the initial Hamiltonian, e.g. H{AiIj) = E{AiIj), V Aip ^ 0. We are 
interested only in those operators A whose action on i^^(R) produces functions in L^(R); 
hence we can choose = cq Atp^ with cq a normalization constant, and we will have 
E = E, \f ip = coAip 7^ 0,'?/' G L^(R). Therefore, the intertwining operator A transforms 
solutions ip of the initial Schrodinger equation into solutions ip of a new Schrodinger 
equation, both equations sharing the initial eigenvalues. However, the new Hamiltonian 
could have a finite number of additional eigenvalues; if this is the case, the corresponding 
eigenfunctions have to be elements of the kernel of A\ as we shall see below. We will 
show now how the intertwining method does work for the simplest case where A is a first 
order differential operator. This case is specially interesting because it provides a general 
framework in which various apparently different methods are included. 
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2.1 The first order intertwining method 

Let the operator involved in the intertwining relationship be the following first order 
differential operator 

A = a= ^ + f3{x). (3) 
ax 

By introducing (|l]) and (|) in (||) we get 

-(3'{x)+f3\x) = V{x)-e, (4) 
V{x) = V{x) +2(3'{x), (5) 

where e is an arbitrary integration constant, and the prime denotes derivative with respect 
to X. Notice the implicit dependence of P{x) and V{x) on e. 

From equations (|l]) and it becomes apparent that H — H = 2/?'(x), from we 
get [a, a = 2(3' (x), therefore (|) and (|) lead to: 

H = a^a + e, H = aa) + e. (6) 

Equation (|^) shows that the first order intertwining relationship (|^), with A given by (0), 
leads in a natural way to the factorization of the intertwined Hamiltonians for each value 
of e G R, provided that [3{x) satisfies with a given V{x). The traditional approach 
to the factorization method involves just a particular solution to the Riccati differential 
equation @) for a specific value of e 0. However, it is well known that a particular 
solution to a Riccati equation allows one to get the corresponding general solution by 



means of two quadratures |jT8[. Mielnik's approach takes this into account [|1^ and gives 
a term in (3{x) unnoticed by the traditional factorization method. With the exception 
of the harmonic oscillator and the hydrogen-like potentials, the general solutions to the 
corresponding Riccati equation remain unexplored for almost all the Infeld and Hull 
cases |1T9|, pO|. An interesting discussion on the relationship between Mielnik's method 



and Darboux transformation was given by Andrianov et al in 1985 |jTT| (see also [20]|). On 



the other hand, equations (^) and (j^) are basic in the construction of the SUSY partner 
Hamiltonians = H — e = a^a, and = H — e = aa^ in supersymmetric quantum 
mechanics. In this language, it is said that equation (|^) relates the two SUSY partner 
potentials V{x) and V{x), and the function j3{x) is named the SUSY potential 0. Notice 
that the two Hamiltonians are positive semi-definite, i.e. (H^) > 0; hence we have 
E > e and E > e. 

If the eigenfunctions and eigenvalues {E} of H are known, the corresponding 
eigenfunctions of H are given hj = {E — e)~^/^ atp with eigenvalues E = E. They form 
an orthonormal set 



iij',^) = [{E' -e){E- e)]-'/'{aij',a^) = = 0, + ^. 

However, they do not span yet the whole of i^^(R) because the solution to (^^l)^ = is 
orthogonal to all the members of the set {'?/'}, {il^eii^) = {E — e)~^^'^{a^-ipe,4') = 0. That 
first order linear differential equation is immediately solved: 

cxe/'^(")''". (7) 
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From it is easy to see that Hip^ = eip^- Hence, ip^ is an eigenfunction of H with 
eigenvalue e. In the case when ip^ is a square integrable function [i.e. ip^ has physical 
sense as a wavefunction) , it has to be added to the set {^p = {E — e)~^^'^a'ijj} in order to 
complete a new basis on L^(R). The set of eigenvalues of H is given by {E} = {e, E}, so 
that the Hamiltonian H has the same spectrum as H plus a new level at e. Indeed, since 
the general solution (3{x) of the Riccati equation has a ^free parameter' A, labeling 
also V{x) (see equation (||)), the Mielnik's method allows one to construct, for a specific 
value of e, a one-parametric family of potentials V{x) which in general will be almost 
isospectral to V{x). 

We are going to disscus now how the Mielnik's method can be further generalized. 
First, let us remark that the factorization energy considered by Infeld-HuU and Mielnik 
is not the only posibility leading to a solvable Riccati equation (^). The key point arises 
when the non-linear first order differential equation (^) is replaced into a homogeneous 
linear second order differential equation by means of the transformation|^: 

P(x) = -^lnu{x), (8) 
which leads to the second order differential equation: 

-^u{x) - [V{x) - e]u{x) = 0. (9) 

If the factorization energy e belongs to {E}, then the solutions u{x) to (|^) are the same 
as the previously known eigenfunctions iplx) of the initial Hamiltonian H. Following 0, 
we shall consider just the cases where the factorization energy e does not belong to the 
spectrum of H. Therefore, the solutions to (|^) do not have direct physical meaning, but 
we have just seen that they naturally lead to the factorization of the Hamiltonians H 
and H in the spirit of Mielnik's approach, providing as well the explicit form for the new 
potentials V{x). 

Let us remark that n different factorization energies ei, €2, e„, lead to n different 
Riccati equations @ and to n one-parametric families of potentials Vi{x), i = 1,2, ...,n, 
almost isospectral to V{x). It must be clear that there will be n different factorizations 
of the initial Hamiltonian H (see equation (^)). With the aim of iterating later the first 
order intertwining technique let us introduce a slightly different notation. Let us write 
e = e'-^'', where k labels the specific real number e ^ {E}, which we will use as factorization 
energy. Moreover, we shall write explicitly the dependence of P{x) and V{x) on e^*^^ (see 
equations (^) and (|])) by making P{x) = P^'^^x) and V{x) = V^''\x). For the sake 
of simplicity the other functions and operators labeled by a tilde will be represented as 
H = H^''\ E = E'^''\ etc. Finally, the number of indices will denote the order of the 
intertwining transformation (Q), as well as the number of '/ree parameters^ A which will 
label the intertwined potential. For example, if'-'^-* represents a first order intertwined 
Hamiltonian whose corresponding one-parametric family of potentials V^'^^x) is almost 
isospectral to V(x), while represents a second order intertwined Hamiltonian with 

^The inverse transformation of (||) is u{x) = exp(— / l3{x)dx). 
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the corresponding two-parametnc family of potentials V^''^\x) almost isospectral to V{x). 
The benefits of such notation will be clear in the next sections. 



2.2 Second order intertwining method 

Let us consider the second order differential operator 

B('=-) = ^+r,(x)^ + ^(^), (10) 

where r]{x) and 7(x) depend implicitly on e^'^-' and e^™^ Introducing if^'^™) = —cP/dx^ + 
yikm)^^^ and (|I0|) in the second order intertwining relationship j^ikm) _ j^{km)j^^ 

with H given in (|l]), and by ordering at different powers of the operator d/dx, we have 
for the corresponding coefficients 

V^'''^\x) = V{x) + 27]'{x), (11) 

2-f{x)=ri^{x)-ri\x)-d-2V{x), (12) 

V"{x) + r]{x)V'{x) = 27(x)r/'(x) - 7"(x), (13) 

where d is an integration constant. Following [|12], we can express the functions V{x), 
^(fcm)(^^ and 7(a;) in terms of r]{x). The substitution of (^) into (0) produces a third 
order differential equation which, after multiplying by ri{x), becomes reduced (by a first 
integration) to the second order differential equation 

VV" - + (^27 - r/' - ^ ^ ^ 2c = 0, (14) 

where we have used once again the equation (|12]) and c is a new integration constant. 
In order to solve the non-linear second order differential equation (p!^, let us try the 



ansatz 21 



T]'{x) = ri'^{x) + z{x)ri{x) + b, (15) 

where the constant b and the function z{x) are to be determined. Introducing ( plSj) in (|1^ 
we get b± = ±2y/c, c > 0, besides two Riccati differential equations 

z'^{x) + ^^~2V{x)-d-b± = 0. (16) 
The identification z^{x) = — 2/5fc(x), z_(x) = —2Pm{x), and 

= c= r , e(-)^6(^), (17) 

allows one to write the solutions of in the form 



where P^"^\x) and I3^^\x) are both solutions to (|^) with factorization energies e'-'"-' and e*^'^-' 
respectively. Equations (|18D and (0) determine the second order intertwining operator 
Notice that, by means of this process, we can derive solutions to the non-linear 
second order differential equation ([I^ ) by just solving the easier non-linear first order 
differential equation (|^) for two different values of the factorization energy e. 

The action of B^^"^^ on the eigenfunctions {ip} of H gives the basic set of eigenfunc- 
tions oc of the intertwined Hamiltonian if^*^™-). On the other hand, the 

solutions to the second order differential equation S^*^™) v]/^™'' = Q are orthogonal to all 
the ^C^*") oc B^'^'^^ip. Now, because the kernel of B'^^^^'^ is a 2-dimensional subspace, the 
maximum number of linearly independent elements is two. Therefore, we are looking for 
two square integrable linearly independent solutions to i^(*:™)^vi/^™) = q such that they 
are simultaneously eigenfunctions of H^^^\ with eigenvalues to be determined. 

Let us now write ^2""^ = Co exp(/ as a generic kernel element of B^ > ; 
then the equation ^^^''^^^^J™^ = can be rewritten as 

^f{km)r. _ ^f^^f{km)r. ^ (f^^\x)f - 7]' (x) + ^(x) = 0. (19) 

ax 

Using equations (^), (|1^) and (pTSf), the general solutions to (|19|) can be obtained: 

<f^^"'\x) = Cu ,..; \ + C„^--rJ^, 20 

W{k,m) W{k,m) 

where u'^^\x) and u^'^^x) are the two t(nj)/i?/5zca/ solutions to (H) with eigenvalues e'^'^^ and 
g(m,)^ respectively. The function W{k^m) = u^''\u^"^^)' — (u^'^^Yu^'^^ is the wronskian of 
u^''\x) and ■u*^™)(x). 

By using (H) and (||) we now rewrite ([l^) as 



/ N ^ (m) (k). u'-''\x)u'-"'\x) d , , 



Notice that the right hand side of (|2T| ) corresponds to the Crum's determinant [17|. 
Hence, the coefficient of the second term of B'^^^^ in ([To|) can be constructed either by 
using (|2T|) or simply by ([T8|). Using now (|T6|), (H) and (|), with rj{x) as given above, it 
is easy to show that the function u'^^^ /W{k^m) is eigenfunction of if*^'^™) with eigenvalue 
^(m)_ similar procedure shows that u^'^'> /W{k^ m) is also an eigenfunction of ifC^"^) with 
eigenvalue e*^'^-*. Taking this into account we write 

where the subindex (cfe) indicates the corresponding eigenvalue e*^™"-* (e*^*^^). Therefore, 
the eigenfunctions of H^^"^^ are {^^^"'\ ^g""^} U {^C^'") oc S(*^™V I B^^'^^ ^ 0}, and the 
corresponding eigenvalues are {ii^C^™)} = {e^'^\ e^^\ E} . Thus, H^^"^^ has the same spec- 
trum as H plus two new levels at e^*"-* and e'-'^-*, provided that ^^^J^™^ and ^1'^^"^^ are square 
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integrable functions. As (|T8D involves two general solutions to the Riccati equation (^), 
thus (x) represents a set of iwo-parametric families of potentials almost isospectral 
to V{x). 

Notice that we have not placed any ordering to the levels e*-"^^ and e'-'^^ so that it 
is possible either that e^"*^ > e'-'^-' or that e*-™"-* < e^^-*. This is implicit in the fact that 
ri{x) in ([T8|) is invariant under the change k ^ m. This will be important in the next 
section, where we are going to show that the results found by means of the second order 
intertwining technique can be found also by the iteration of two successive first order 
intertwining transformations. 

2.3 Iterative factorization 

Let us come back to the results of Section 2.1 and suppose that by means of the first 
order intertwining technique we have obtained a Hamiltonian if^^^ from H using a certain 
factorization energy e'-'^^ We look now for a new Hamiltonian /i^'^™) = —cP/dx^ + v'^^^\x) 
attainable from H^^^ by means of the first order intertwining method, i.e., h^^^^f and if'-'^^ 
satisfy the following relationship 

where a*^'^™^ is the first order differential operator a'^'^™) = d/dx + As usual, we 

arrive at the standard equations linking V^^\x) and P^^"^\x): 

- + = -em = 2(/3('=))' + V{x) - (23) 

= V^^\x) + 2(/?('=™))'(x) = V{x) + 2[(3^'''^\x) + /^(^Hx)]', (24) 

where Sm is a new integration constant (factorization energy) and we have used equation 
d^). The intertwined Hamiltonians once again become factorized 

^(fc) = ^(fcm)t^(fc„^) ^ j^^km) ^ ^(fcm)^(/cm)t ^ 

The eigenvalues of /i^^™) are given by = {e:^, and the eigenfunctions by 

^^fcm)y{^(fcm) ^ ^^(k) _^^y 1/2 ^(km) ^(k) | ^(fcm)^(fc) ^ Q}, where the missiug statc ^^^'"^ cx 

exp(/ I3^^'^\x)dx) corresponds to the eigenvalue Sm- 

The key point here becomes to find the general solution to the new Riccati's differential 
equation (^) for some Em G R. The obvious solution (/3('^™))'(x) = — (/3'^'^^)'(x), for which 
Em = corresponds to a Hamiltonian h^''"^'> equal to H because v^''™'\x) = V{x). 
Avoiding this trivial solution we shall consider Em 7^ Non-trivial solutions to (p3D 
can be found by inspecting the following intertwining relationship: 

hikm)^{km)^{k) ^ a^km)jj{k)^{k) ^ a'^km) ^{k) ^ ^25) 

Equation ( P^D shows that /i^^™-) and H are related by the product of two first order 
intertwining operators, i.e., by the second order differential operator C^^"^"^ = 
An interesting question immediately arises: is the Hamiltonian h^'^™'^ the same as the 
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jj(km) Section 2.2? If = B^''"'^ the answer is affirmative, and in this case the 

SUSY potential P^''\x), the iterated new SUSY potential P^''"^\x), and the functions 7(x) 
and r]{x) are related by the equations 

r^(x) = (3('''^\x) + f3'^''\x), 7(x) = -^(3^''"'\x) + (3^''"'\x) (3^''\x) . (26) 

dx 



Comparing (0) and (|26|) we get the solutions of (^3]) (with Er 



r(ni)' 



^""•'W--^'"W-i /^.^>W:;3.>.M j- 

where /?(*^)(x) (/^("^^(x)) is a solution of (^ associated to the factorization energy e'-'^-' (e'-"^-'). 
Equation ( p7|) is nothing but the theorem of Fernandez et al The eigenf unctions of 
j^{km) _ jj{km) ^^^^ given by the basic set 

^(fcm) ^ ^^{k) _ ^(m)^|-l/2^{fcm)^(fe)^ ^(km) = gM)-l/2(^_g _ g W ) " V2Q(fc"^)a 

plus the new missing state 

Cm ' 

which, by definition, is such that a^''"^^'^'^!^''^^ = 0. Hence, by construction _B('='")t\[/('«™) = 
Q(fc)tQ(fcm)tv|/(fcm) _ g_ ^Yie other hand, remembering that a^'^^^tp^'^^ = 0, we get 

Now, by using (^6]) , (H) and after (|2l|) , we can write 

which corresponds to the definition of \l/(^™)(x) given in (^21). A similar procedure but 
now using (pT]), (||) and (|21|), allows one to recover the corresponding expression for 
^(fcm)^^)_ Finally, the set of eigenvalues of /i^^™) = i/^'^™) is given by {^(*^™) = E^'^™)} = 
{e('^),e(^),E}, just as we have shown in Section 2.2. 

It is worth to notice that the operator ijC^™) = a^'^^^a^'^'-' does not factorize the inter- 
twined Hamiltonians H and H'^^^\ but some of their quadratic forms: 

= (iJ-e('=))(/7-e('")) = (if-e(™))(if-e(^')), (28) 

and 

j^{km) j^{km)\ _ ^jj(km) _ ^(m)'j^jj(km) _ ^ik)y ^2Q'j 

Up to now, from the general solution to (^) for each factorization energy e^''\ we 
have derived a one-parametric family of potentials V^''\x) isospectral to V{x) but by a 
new level at e*^'^-'. From two of these solutions with e*^'^-' ^ e^"*^ we have derived as well 
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iwo-parametric families isospectral to V{x) but by two new levels at e*^™-* and e^'^-'. It is 
possible to iterate further the first order method by considering now the general solutions 
to (^), associated to three different factorization energies e'-'^-*, e^"^\ e^^-*. The key Riccati 
equation to be solved becomes now: 

The algorithm ( p7|) allows one to write (3'^^'^^\x) in terms of two different solutions of 
which leads to the following 3-parametric family of potentials V^'-'^™^-* (x) : 



^(fcmO(a;) = + 2(/?(^™'))'(x) = Vix) + 2[/?(^™')(a;) + P^'''^\x) + P^''\x)]', (30) 

having the same spectrum as V{x) but by three new energy levels at e^'), e^""^ and e^'^) [|15[. 
The corresponding third order intertwining operator is given by C'^^'"^) = a'^^'^^^ a'^^^^ a'^'^\ 
This leads to the factorization of some six order differential operators, related to the two 
intertwinned Hamiltonians by means of: 

Q{kml)^ (j{kml) ^ _ gW^^ff _ e(''"))(if — gC^)) 
^{kml) ^{kml)^ ^jj{kml) ^ (kml) ^(m) ^ (kml) g^^)) 

This process can be continued at will, and the corresponding formulae for the key functions 
can be algebraically obtained from the solutions P^''\x) to (H). It is interesting to see now 
how the process does work explicitly for a physically meaningful system. 



3 The hydrogen-like potentials 

Let us consider a single electron in the field produced by a nucleus with charge Ze, where 
Z is the number of protons in the nucleus^. The time-independent Schrodinger equation 
is given by: 

_ V2^(f) - -^(r) =E^(f), (31) 
r 

where we are using the units of = fi^/Ze^m, and £ = Z/2rB for the coordinate 
and energy respectively. Due to the spherical symmetry, the standard choice \l/(r^ = 
R{r)Y{6, (f)) separates ( PT| ) into a radial and an angular equation. The angular solutions 
are given by the spherical harmonics Y{^{9,(j)) and the energy spectrum can be finally 
obtained by solving the radial equation: 



cP lil 
+ — 



dr 



ip{r) = Eip{r) 



(32) 



= 1 describes the Hydrogen atom, while ions or atoms with only one electron in the outermost 
shell can be considered as hydrogen- like systems with Z > I. 
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where, / = 0, 1, 2, is the azimutal quantum number and, for the sake of simphcity, we 
will work with ^/'(r) = rR{r), < r < +oo. The set of solutions of (|32D span the whole 
of L^(R~^), with an inner product defined by {ip,'^'') = 47r Jf^°° ip{r)ip'{r)dr < oo. The 
differential operator of the left hand side of (^) will be refered as the radial Hamiltonian, 
and it will be denoted by Hi, where 

m = -^^+mx); ^(x)-^-f (33) 

The eigenvalues, for a fixed /, are given by the well known formula 

1 

(/ + K) 



En = EiK = -j^-^, K= 1,2,3,... (34) 



3.1 The new one-parametric families 

Let us consider the first order intertwining relationship H^^^a^l''' = a\''^Hi where, as a result 
of the further calculations, we are labeling from the beginning the unknown Hamiltonian 

(k) 

Hi_[ with the subindex / — 1 because it will arise a centrifugal term for the new potential 
with exactly that index. In order to solve the corresponding Riccati equation (^) we shall 
use its equivalent second order differential equation (P) rewritten as 



rf2 ^ /(/ + 1) 2 



«f)(r) = 6pM'=)(r). (35) 



Comparing with (R^) and (B^, it is natural to propose the factorization energies as: 



(l + k) 



^ ky^K, l>0, (36) 



2' 



where, in general, k could be either a discrete as well as a continuous real number. On 
the other hand, note that the case k = K reproduces the well known solutions to the 
physical eigenproblem ( p^ - |3^ ); hence, following Mielnik we have taken k ^ K. It is clear 



that, for a fixed value of I, we have ep'' ^ Eix, y k ^ K. 
In order to solve ( ^5]) we make the transformation 

«p)(r) =r-'e'-/('+'^)<l>p)(r), (37) 

leading to a confluent hypergeometric equation for $P''(r), whose general solution for the 
discrete values k = 0,-1, —2, ...,—(/ — 1) is given by the linear combination of confluent 
hypergeometric functions [^, |2^ : 

$W(r) = ^Fi[k,-2l,-2r/{l+k)]-iyik[2r/{l+k)]^+^\F,[l+k+2l,2+2l,-2r/{l+k)], (38) 
with 

r(l + |A;|) , , , 

- r(2/ + 2) (-2/),| - ' > »• '^«) 
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(k) 

where A; is a constant to be determined and 



;-20ifc 



r(-2/+ 



[-21 + \k\ - l)(-2/+ -2) 



Note that $ 



r(-2/) 

0) = 1, while its asymptotic behaviour is given by 



-21). 



2r 



\k\ ^ _ ^k) 



I - \k\ 



:-2i) 



\k\ 



Hence, the solutions ( P7| ) are divergent at the origin as well as in the limit r — »• cxo, and 
it is clear that u\''\r) ^ L^(R+). This is the reason why we consider them as unphysical 
solutions to the Schrodinger equation (|3^) with the atypical eigenvalues (|36D . Let us 
remark that the discreteness of k in ( |57| - PUD leads to the most transparent case available 
with the mathematical tools at hand [Q. The general solution to the corresponding 
Riccati equation (H) arises after introducing ( p7|) in (H), which gives: 



I 



1 



Using (H) and 



r (l + k) 

we get the first order intertwined potentials 



/ > 0. 



(40) 



Vl':i{r) = Vl^^ir)-2—\n^r>{r), 



I > 0. 



(41) 



The second term of (|4l|) is free of singularities if 

[—oo, 1), for \k\ even; 



A 



(fe) 



(42) 



;i,oo) 



for I A; I odd. 



Hence, the new potentials Vi^^l{r) have the same singularity at r = as V5_i(r), 

provided that A^"* takes only values in the domain (^2]) . Since the second term of Vi^_!^l{r) 

tends to zero when r ^ and r — oo, we conclude that V^^i(r) behaves as V5_i(r) at the 
ends of [0, oo). 



Now, the eigenfunctions of are obtained from the action of a\'^' on the eigenfunc- 
tions {ipni{r)} of Hi, ipnl-i = i^n — O'T'^ i^niir) , plus the isolated giomid state: 

j,l ^-r / {l+k) 



(k) 



where $P''(r) is given by 



Ik 



and 



Ik- 



$f)(r) 



(43) 



;-2/)|fe| / (2/)! 



1/2 



-1,-2,..., 



-(/-I). (44) 
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It is important to notice that, despites the resemblance between the missing states ( ^3[ ) 
and the unphysical functions (|37D, "ipi-i^^k ~ V'^P^ J^st the first ones become square 
integrable functions. The eigenvalues of -ff/^i (for fixed / and k) are: 



By comparing with (|3^ ) it is apparent that our families of potentials ( |41| ) have, in general, 
the same spectrum as the corresponding spectrum of the hydrogen-like potentials V;_i(r), 
plus a new level at e^^^^ 

3.2 Particular one-parametric families 

Notice that, for fixed /, k can take / different values k = 0, —1, —2, —(/ — I). Hence, there 
are / different factorization energies e[^^ generating / non-equivalent families of solvable 
potentials V^!^i^(r). In particular, let us notice that e;^'' = = L{1), with A['^'* = 0, 

lead to the standard factorization of the hydrogen-like potentials as presented by Infeld- 
HuU (see equation 8.1.3, and page 68 of [@]). Thus, our technique can be considered as 
a generalization of the Infeld-HuU factorization where, instead of looking for particular 
solutions to (^ with ef'^ = L{1), we look for the general solutions with ep = L{1 + k) . In 
a similar way, our results can be seen as a generalization of those derived by Fernandez, 
who found a general solution to (^ but just in the case with = Hence, for a fixed 
/ and taking k = 0, the equation ( ^Sj) leads to E^^^ ^ = Ei_i^k, ^-c, Vjf°| (r) is a family of 



potentials strictly isospectral to Vj_i(r). On the other hand, for k ^ 0, we have not only 
e[^i\ 7^ Ei_i^K, but e[^i\ < Ei_i k=i = + 1)^- In this case the spectrum of Hf^'^'^ 

is almost the same as the one of the difference resting in the ground state energy 

level. It is clear now that the functions (^) represent a set of / — 1 families of potentials 
almost isospectral to V/__i(r), plus a family strictly isospectral to V5_i(r). 

3.2.1 The case A; = 

The member of the family (|41|), for /c = 0, is given by 

^SM = ^^-iW-2^1n$r(r), />0, (46) 



where 



2, 2i+l >(0) 



From (|i2| ) we see that Vi^\{r) does not have more singularities than V5_i(r) if A^^^* G 
(—00,1). In particular, for A['^'* = 0, we have <^f\r) = 1 and V5*^](r) = V/_i(r), which 



means that V;_i(r) itself is a member of the family (|46|) . On the other hand, we get also the 
standard SUSY potential W{r) = j3f'\r) = l/r — l/l, besides the corresponding first order 
operator af'^ = d/dr+l/r — l/l, both of them used to solve the hydrogen- like Hamiltonians 
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Il], H, H]. Let us remark once again that V^L°J(r) with X\ 



(20!(V2)^'+Wr^ leads to the 



one-parametric family of isospectral hydrogen-hke potentials derived by Fernandez IQ. If 
/ = 1 and Xi'^ 1, the potential (^6]) tends also to the one found by Abraham and Moses 



3.2.2 The case k = -1 

Let us take now k = —1. Thus, the potential (|41|) becomes 



where 



V, 



Z-1 



1(1-1) 



(-1), 



r), Z>1, 



1 + 



A 



(-1) 



(2/-1)! V/-1 



21-1 



^2l^-2x/(l-l) 



[/(/ 



dx 



x\ 



(48) 



(49) 



In this case A[ G (l,oo), and although V;_i(r) governs the behaviour at the origin 
and at infinity of Vi^_i\r), it is not a member of the family (P8|). The corresponding 
eigenvalues are given by 



(-1) 

l-l,K 



(/_1)2' (/ + 1)2' (/ + 2)2^ 



/ > 1. 



(50) 



Notice the 'jump' between the first two energy levels of h\_]^] this produces a hole in the 
spectrum ( |50D because the level at is absent. Comparing (|50|) with (|3^ it turns 



out that ViO^\r) represents a family of potentials almost isospectral to V5_i(r). As a 
particular example let us consider 1 = 2: in this case Vi~^\r) has the same levels, — 
as Vi{r) for n > 3. However, the ground state energy level of Vl~^^ (r) is at —1, which 



is forbidden for Vi(r) whose ground state is at —1/4. The potentials ( ^^ - |49| ) have been 
recently reported by the author (see equations (10) and (22) of fl^) and they are plotted 
in Figure 1 for / = 2 and different values of A2~^''. 



3.2.3 The case \k\ > 1 



(k) 

For \k\ > 1, there is a gap between the two first energy levels of Hi_[, which grows up as 
I A; I increases. For instance, when = 2, and / = 3, the levels labeled by n = 2 and n = 3 

(—2) 

are not present in the spectrum oi H2 .A similar result, but with a different method 



and Hamiltonian, has been derived by Samsonov [11]. 

Recently, the SUSY potential uj{x) = a/7 — 7/x, with a > and 7 > 0, has been used 
to construct a 'conditionally exactly solvable (CES) potential' V^{x) for the Hydrogen 
atom problem p5|, |2^. In such a paper a is a coupling constant while 7 is a parameter 
taking some strange values as 7 = 2.8 (see Fig. 7 in p6i)- Although there is a close 



relationship between that potential V+{x) and our first order intertwined potential VJ 
it is important some caution: 
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z) If 7 > is not an integer, the SUSY potential u(x) of I^G] is not related to the 
physically relevant SUSY potential W{r) for the Hydrogen atom, and in this case 
V+{x) is the 'CES potential' for some mathematical problem. If the so called CES 
potentials derived in are to be related with a "physical potential" , the condition 



7 > is not sufficient. 

ii) If 7 is an integer greater than zero, the SUSY potential u{x) can be related to 



the Hydrogen atom problem. In this case the function (5.14) of p6| has a term 
iFi(— 7 — a/p, —27, 2px), which is meaningless at 2px = unless 7 + a/p is an 
integer greater than or equal to zero [|^. Indeed, since 7 is a positive integer, the 
fraction a/p is an integer such that 7 > —a/ p. Hence, the condition a > is neither 
sufficient to ensure the generality of the potential V+{x). 

Taking this into account, a general 'CES potential' V^{x) really related to the Hydro- 
gen atom problem would take for instance —a/p = 7 + /c, where k is an integer such that 



A; < 0. In particular, for 7 = /, and x = r/a, the CES potential V+{x) derived in p6 



would be equal to our first order intertwined potential Vi^^l{r), provided that uik = (3/a, 
where a and /3 are constants. 

We are going to derive next some two-parametric families by means of the technique 
discussed in sections 2.2 and 2.3 for the particular case we are dealing with. 

3.3 The new two-parametric families 

In Section 2.2 we have seen that the second order intertwining relationship H^^^ b\^^ = 
Bi'^^^Hi provides new solvable potentials if we are able to find some solutions ri{x) to 
equation ([Tl|). Then we also have shown that the general form of that solutions is given 
by: 

/ \ / W M\ ) ^ ^ ( ^PH'") 1 1 (fe) / (m) 



dr 

\ "'I V 

For our particular system we take the functions uf^\r) as given by (P?!). The behaviour 
of our rjir) at the ends of [0, 00) is given by 



2l + k + m (1-2/) 
{l + k){l + my 



Hence, it is natural to write the new potential Vi^'!^\r) as (see equation (|Tl|) ): 

v}^:^\r) = Vi.,{r)+2a\r), (52) 

where a{r) = r]{r) + (1 — 2/)/r is an appropriate function making evident the limit 
Vi^_!^^\r) — s> Vi_2(^) when r +00, or r ^ 0. The two-parametric domain of A['^'* and 
Xf^\ for which a'{r) is free of singularities, is determined by the parities of k and m. In 
order to make transparent the choice of that domain we take the following convention: n 



14 



is the index labelling the factorization energy e^""^ defined by tf"' = ma.x{e\''' , e^'), and s 



is such that e\^^ = min(e['^'*, e|™^). Hence, we have: 



A;""* G (—00,1), A;"*"* G (— cxD, 1) ; \n\ even, |s| odd 
A;""* G (1,00), A;'*'' G (1,00) ; \n\ odd, \s\ even. 



(53) 



\i G (— cxD, 1), A^"* G (l,C)o) ; \n\ even, |s| even 



Ar G(i,oo), Ar G(-oo,i 

The eigenfunctions of h\^^ are given by 



\n\ odd, Isl odd 



(54) 



(i?n-6S'"))(i?„-el 



(55) 



and 



_{k) (m) "/ 



(m). 



(fcm) / \ p / [I') 



(56) 



where Cik is given by (|^). The corresponding eigenvalues are: 



Tp(km) 
^l-2,K 



m] 



(/ + ky 



U 



(l + Ky 



K>0. 



(57) 



Let us remark that the only restriction for the factorization energies, ej 7^ i 
leads to the fact that ri{r) and V5!^;^'*(r) are symmetric under the change k —>■ m and 

However, this symmetry is broken for the intermediate 
jg obtained after two iterated first order 



T r(km) / \ T r(mk) / 

Vice versa, Vi_2 [r) = Vi_2 v 
one-parametric potentials arising when VJ_2 
intertwining procedures. In order to see that, let us suppose that e 
two different ways to get the potentials 



(fc) 



> e 
(fc) 



(m) 



We have 



we can add first the level e] , and then the 
level e[™'''; the other option is to perform the two operations in the opposite order. 

(I) In the first case, after the first step we shall get the Hamiltonian Hi^^^, with 

(k) 

a spectrum equal to the spectrum of Hi plus a new level at e[ < Ei^k=i, where 
Ei^K=i = — + 1)^ is the ground state energy level of Hi. After the second step. 



the final Hamiltonian Hf^'' has the same spectrum as -ff/^i plus a new level at 



{k) 



e[™^ < e^^ ■ In both of these procedures the factorization energy (the ground state 
energy of the new Hamiltonian) is less than the ground state energy of the initial 
Hamiltonian, just as it happens in the SUSY approach. 
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(II) Taking now the two operations in the opposite order, the first intertwined 
Hamiltonian Hi^^ has a spectrum equal to the one of Hi plus a new level at 
^ Ei^K=ii in agreement once again with the SUSY approach. However, the 
next intertwined Hamiltonian h1^2^ is derived by adding a new level at ef^'^ > e\'^\ 

i.e., above the ground state energy level of h\^1, which is against the SUSY doctrine 
but it is valid in our treatment. 

A final point is that the equality of the two Hamiltonians -ff/^™'' and Hi^2 \ gotten of 
two iterations of the first order transformations, leads to two different factorizations of 
the global second order intertwining operator 

n(fcm) _ (km) (k) _ (mk) (m) _ n(mk) 

3.4 Particular two- parametric families 

The simplest two-parametric family of potentials (^) corresponds to A; = — 1 and m = 0. 
In this case, by means of ( |5lD we get the solution to (p^: 

1-2/ i d ^ f e'/Ki-i)^(-^)fr) , , , 
^^^^-W^^[d-r'^[ ^(°)(.) ' 

which has been recently reported in ||1^, where ^f'\r), and $[ ^^(r) are given by (^Tf ) and 
(|49|) , respectively. The appropriate tt(7o-parametric domain (see (|53D and (H)) is given by 
A["^\ A{°^ G (-00, 1). From and (0) it is clear that 1^^'°^ (r) is strictly isospectral 
to V^_2(r). Notice that a{"^^ = A{°^ = leads to ■^L~2^'°V) = V/_2(r). Moreover, for 
Af^ = and X\~^^ = [{I - l)/2f-\2l - 2)!/7;_i, the potentials (H) become the family 
( |i6| ) derived by Fernandez, provided I is changed by / + 1. We make concrete now the 
discussion of the cases (I) and (II) of Section 3.3. 

(I) Departing from the Hamiltonian Hi, by means of a first order intertwining op- 
erator we get the Hamiltonian Hi^^, whose potential is given by (^61), where we 
have added the new level ef'^ = below the ground state energy level of Hi, 

Ei^K=i = — 1/(^ + 1)^- The iteration of this procedure leads to h\^2 ^\ whose po- 
tential is given by ([5^ ) and (|5^). In this second step we have added the new level 
e\ = — 1/(/ — 1)^ below the previous ground state energy level ef'^ = —1//^, which 
agrees with the usual SUSY asumption. 

(II) We depart now from Hi to get h\_^ , whose potential is given by (^8]). In 
this step, we have added the new level e^^^"* = — 1/(/ — 1)^ below the ground state 
energy level of Hi but we have left a hole to be filled during the second step (see 
Section 3.2.2). The next step provides now h\Z2'^\ whose potential is given again 
by (13) and (|S|). Notice that = A-T^V) + - A~^\r)- The 
new energy level ef^ = -l/f is now added above the ground state energy level 
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e| = — 1/(/ — 1)^ of Hi_^ but below the ground state energy level + 1)^ 

of Hi. This last procedure fills the hole generated during the first step of the whole 
procedure which is not typical in the usual SUSY asumption. 

A different choice of k and m produces families whose spectra are almost equal to the 
corresponding hydrogen-like spectra, with the two first energy levels different from the 
levels of the corresponding hydrogen-like potential. For instance, ii k = —3, m = 0, and 
/ = 4, the potential V2 ^ (r) has a spectrum of the form — but without the levels 
at n = 2 and n = 3, similarly as it happens in the first order case mentioned in Section 
3.2.3. This is illustrated in Figure 2. A more exotic case is obtained if = —4, m = — 1, 
and / = 5, where the levels with n = 2, n = 3, and n = 5 are not present for the potential 



4 Concluding remarks 

The nth-order intertwining technique developed in this paper is intended to derive n- 
parametric families of potentials which can be isospectral, or almost isospectral to some 
well known potential. In particular, we have successfully derived new potentials from the 
hydrogen-like atom for the one parametric and two parametric cases. We have shown also 
that the derivation of multi-parametric families of almost isospectral potentials becomes 
mathematically easy by the iteration of the first order intertwining technique. Although 
the only trouble in our approach could be to find the solutions to a set of n Riccati type 
equations, characterized by n different factorization energies, we have shown that the 
eigenvalues problem for the initial potential allows one to get that solutions by means of 
unphysical eigenfunctions. Our one-parametric results for the hydrogen-like atom show 
how the technique effectively works, moreover, our resulting families of potentials are 
more general than other results previously derived |^, ^, On the other hand, our 
two-parametric potentials represent a further generalization of the one-parametric cases 
presented in this paper. Finally, we have shown that the construction of a SUSY partner 
Hamiltonian using a factorization energy e less than the ground-state energy of the depar- 
ture Hamiltonian is not a general rule in SUSY QM. In particular, the generation of holes 
in the spectra of the n-parametric families of almost isospectral hydrogen-like potentials 
becomes the key-stone in the construction of the atypical n-SUSY partners of the Hydro- 
gen potential derived in this paper. These techniques can be applied to other physically 
interesting situations, e.g., for systems with continuous spectra as the free particle. A 



detailed study in this direction will be given elsewhere |27|. 
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1 2 3 4 5 

r 



Figure 1: Some members of the one-parametric family of potentials P8| ) for I = 2 and 
the indicated values of A2 ■ The broken curve represents the hydrogen-like potential 
V5_i(r). Notice the arising of global minimum for Vi~^\r), which produces a well which 
can be modified by changing the values of the parameter A2~^''. 
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Figure 2: Diagram showing the first energy levels of H2 and H2 ' ■ Notice the gap 
between the ground state energy level 64"^^ = — 1 and the first excited state 64*^^ = — 1/42 
of the tt^JO-parametric family of Hamiltonians H2 . 
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